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ABSTRACT 

Simulation generated probabilities for an Anti-Submarine 
Barrier composed of submarines are examined from a statisti- 
cal viewpoint. A probabilistic model which is not generally 
considered to be applicable to this case is demonstrated 
to be statistically supported. The applicability of the 
model is justified probabilistically. A statistical esti- 
mating relationship is then developed to estimate the sole 
input parameter from submarine FIGURE-OP- MERIT. 
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I. INTRODUCTION 



A. THE TACTICAL PROBLEM 

The encounter treated is that of an Anti-Submarine 
Warfare barrier consisting of submarines. Each barrier 
submarine Is confined to an operating area, hereafter called 
a cell, due to the requirement to remain In a safe haven. 

The cell Is part of a larger barrier across a restricted 
passage of width D. The transiting submarine Is obliged 
to pass through the passage, and hence the barrier. The 
barrier has uniform cells and may be composed of several 
cells In depth. All detections are made with passive SONAR. 

No external Intelligence of the position of either opposing 
.force Is allowed. All detections are made by the submarines 
and all kills are made without the aid of external forces . 
Transiting submarines are alone, and the time between tran- 
sits Is sufficiently long to prevent there being any effect 
of one transit upon another. If cells are arranged In se- 
quence along the transit track, they are sufficiently sepa- 
rated to prevent there being any effect of one engagement 
upon another (This assumption Is relaxed for purposes of 
investigation in Section II.C). The barrier is sufficiently 
wide (D Is sufficiently large) to make any special effects 
from the cells on the ends negligible. 

By virtue of the above assumptions, the encounter may 
be studied as a single cell with a single transiting submarine. 
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and a single barrier submarine. The barrier In total may 
then be constructed probabilistically from the building 
block of one cell. Figure (1) Is a geographical diagram 
of one cell In the barrier, with pertinent variables 
labelled. 

B. THE SIMULATION 

For the purposes of simulation, the following assump- 
tions were deemed appropriate : 

(I) The transiting submarine chooses a penetration 
point uniformly along the cell front, since the limits of 
the cell are unknown to It. 

(II) The transiting submarine proceeds at transit speed 
(V,p) until detecting the barrier submarine, at which time 

It slows to quietest speed to evade. 

(III) The transiting submarine penetrates the cell 
front perpendicularly, and continues on a straight track 
throughout . 

(iv) No counter measures are employed. 

(v) The barrier submarine proceeds at patrol speed (Vg) 
across Its cell until a detection, at which time It commences 
interception tactics. For a short Initial period It steers 
the target’s bearing, after which time sufficient Informa- 
tion is available for a firm track. When firm track Is 
established, the barrier submarine proceeds to the Intercept 
point . 
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Figure 1 

The Tactical Scenario 
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The simulation was by digital computer. It serves the 
purpose of providing data for analysis only, and is des- 
cribed here to give insight into the likely characteristics 
of the data. The validity of the simulation is not a subject 

i 

of direct analysis. The ability of the data to support 
reasonable models, and Internal consistency were the facets 
of validity tested. The simulation was constructed and 
run by the Center for Naval Analyses and should be con- 
sidered separate from the analysis contained herein. 

The simulation is a time-stepped, event-store, Monte- 
Carlo program. Each stochastic process is modeled stochas- 
tically, and each probabilistic event is determined by ran- 
dom number. Input parameters Include acoustic character- 
istics of both submarines, sea conditions, initial veloci- 
ties of both submarines and weapons load of both submarines. 
Sea conditions were not varied for any simulation runs. 

Two hundred and fifty replications of each submarine 
pairing were conducted. Taking any particular binary value 
output by the simulation to be a bernoulli random variable 
(l.e.: target killed or not killed), the variance of any 

probability estimator is 

Var (p) = - -P I 

n 

This quantity is maximized for a p of .5, so the maximum 
possible variance of such output is .001. 
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with two hundred and fifty replications the Central 
Limit Theorem will be assumed to hold, and hence the 
estimators will be assumed to be distributed as gausslan 
random variables. 

C. THE ANALYTICAL PROBLEM 

The original analysis of the tactical problem was 
conducted by simulation. Seven distinct barrier submarine 
hull and sonar combinations and four distinct translter 
hull and sonar combinations were Included. 

Various quantities were compiled and computed from the 
simulation run, most for the purpose of checking consistency. 
The sole parameters finally desired from this simulation 
were the probabilities of kill for' the barrier and the 
transiting submarines. 

The number of distinct simulations is expressed as the 
product of (i) the number of cell widths sufficiently close 
to allow interpolation and sufficiently numerous to avoid 
extrapolation; (ii) the number of opposing velocity pairs; 
(ill) the product of the numbers of distinct opposing 
submarines. The essence of the problem is thus to reduce 
the Information contained in such output to a tractable 
analytical expression. The purpose of this reduction is 
four- fold. First, the information must be made convenient 
for a decision maker. Second, the information should be 
expressed as a function of tactical parameters and 
characteristics to allow inspection of the dynamics of the 
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problem and to assess validity. Third, the Information 
should be expressed as a function of Input parameters and 
decision variables to enable the comparison of alternative 
courses of action. Fourth, an acceptable analytical model 
can either reduce or eliminate the need for further 
simulation computer time, should additional range of input 
variables be desired. 
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II. THE SEARCH MODEL 



The following discussion treats the searcher as moving 
and the target as stationary. In order to be in consonance 
with discussions in the literature. Clearly, the inter- 
changeability of searcher and target may be used to achieve 
a desired result. 

A common model in search theory is the Random Search 
Model, now a classic in the field [see ref. 1]. The 
assumptions for'this search law to hold are listed in 
ref. 1 as 

(i) The target (here the searcher) is uniformly 
distributed within A. 

(ii) The observer's (here "the target's") path is 
random in A in the sense of as having its different (not 
too near) portions placed Independently of one another 
in A. 

(lil) On any portion of the path which is small 
relatively to the total length of path but decidedly larger 
than the range of possible detection, the observer always 
detects the target (here "the target is always detected by 
the searcher") within the lateral range W/2 on either side 
of the path and never beyond. This last assumption can be 
relaxed to the extent that a number W (search sweep width) 
may be defined by integrating over the lateral range curve 
P(x) where 
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P(x) = P(detect a target passing 
at lateral range x) 

A searcher with lateral range curve P(x) will detect as 
many targets in the entire plane as one with a definite 
range law detector such as indicated by assumption (iii), 
as long as the targets are distributed uniformly, and as 
long as 

00 

/ P(x) dx = W 
—00 

where W is the diameter of the definite search law detector. 
Under the Random Search Law, 

P(detection) = 1 - 

‘ W = sweep width 

I = track length 
A = area containing target 

If the searcher transits through A then 



£ _ 1 
A “ C 

C = width of A 



thus 



P(detection) = 1 



-W/C 

e 
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For the purposes of this analysis, the assumptions of 



the Random Search Law hold, with the exception of (11). 

The searcher (target In this context) Is not searching 
along a random track, but Is actually proceeding directly 
along an effectively zig-zag non-overlapping track relative 
to the target. 

A- second model Is developed In ref. 1 which may also 
be applied. This model Is a barrier model. It will be 
referred to here as the Channel Patrol Model. 

Suppose a searcher patrols at speed v across a channel 
of width c. Then one patrol (here a channel crossing) Is 
completed. In ^ hours. The target Is approaching at speed 
u straight down the channel. During one patrol of the 

r\ 

searcher he travels ^ u miles. Relative to the target, 
the searcher Is traveling a zig-zag path at relative 
speed w, where 



Thus the total area within which the target may be 
found In one patrol Is 




(see Figure 2). 



A = (height) (base) = (u ^) c 
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Figure 2 

The Channel Patrol Model 



19 



The total relative area swept by the searcher in one 
patrol is 

B = (velocity) (time) (sweep width) 

= (w) (|) W 



If the target is located within this relative area 
uniformly, then the probability of detection is 



P(det) 




w (-) W 
u (f) C 



W w 
C u 



Noting that 



w 



= + v^ 



= uv^ 



1 + 



(V)2 

u 



then 



P(aet) = = w ic 



if ^ k << 1 then approximately 



P(det ) 




This approximation Involves some reduction in the 
probability of kill, since 



1 
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This may be Justified in that the constrictions of a cell 
or channel will deny the barrier submarine the portion of 
his sensor's range which is beyond the limits of his cell. 
This phenomenon will be called clipping, and will be 
discussed in Appendix B. 

If the probability of kill given detection is p, for 
all targets then 



P(kill) = 1 




k 



k = 1 random search 

model 



k = 



channel patrol 
model 



The quantities p and k are statistically nuisance parameters, 
and may be suppressed by defining a new quantity, kill 
sweep width, and also denoted W, so that 

_ 

(1) P(klll) = 1 - e ^ 

W = kill sweep width 

The advantage of this representation is that p is no 
longer needed, provided that kill sweep width can be 
determined, as will be the case. 
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B. MODEL VERIFICATION 



Rather than questioning the reasonableness of the 
assumptions, or the applicability of the model to the 
tactical problem, an alternative approach seems fruitful. 
Let us suppose the applicability of the model and then 
test the supposition statistically. 

It is natural to suppose that W (sweep width) is a 
function of at least (a) speed of both submarines and 
(b) quietness of both submarines. It is explicitly not 
a function of C (cell width). The obvious test of the 
Random Search Model is to Invert the formula, injecting 
known C and monte-carlo estimates of P(kill) to obtain 
W estimates, and then test the variability of w with 
respect to C. Now, 

W 

- T “ C 

Pk - 1 - e 
there fore 

W = -C logg (1 - Pj^) 

Both a W,p and Wg are found, which are those associated 
with the translter and barrier respectively. 

From considerations of the data, not Included here, it 
is immediately obvious that the variation in W with respect 
to C is very small. This impression was verified statis- 
tically using analysis of variance. It was found that the 
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variations in W were Justifiably attributed to statistical 
fluctuation in the monte-carlo estimation of P(kill), and 
not to any systematic model inconsistencies. The detection/ 
kill model is thus supported statistically (see statistical 
appendix) in the sense that it falls to be contradicted. 

A traditional approach to the simulation data used here 
probably involve the estimation of W from technological 
considerations, and a comparison of this computed value and 
the attendant probabilities to the results of the simulation. 
The approach taken herein is to take the simulation results 
and derive the implied values of W. Validation of the model 
would then involve consistency within the model. The main 
reason for this approach is that only kill probabilities 
are desired, whereas typical estimation of w is for search 
sweep width. Thus inversion of the procedure allows for 
elimination of the necessity of analytically solving the 
klnematical problem of kill, counter-detection and evasion, 
which is the primary source of difficulty in the traditional 
approach. Additionally, should the predicted sweep width be 
sufficiently different from an experimentally derived value, 
then the classical approach offers little recourse except 
to discard the predicted sweep width, the experimental data 
or the model. 
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C. IMPLICATIONS OF THE MODEL 



There are very important tactical considerations which 
might argue for various barrier configurations for 
submarines.' Two building block configurations are here 
examined to determine their probabilistic advisability. 

The first is the "stacked" barrier, with submarines deployed 
in depth. The second is the "packed" barrier with 
submarines deployed along a line (see Pig. 3). 

Now, suppose that there is some larger barrier of 
width D in which these cells are placed. Alternative I is 
to place one packed barrier, alternative II is to stack k 
less dense barriers with the same number of submarines, 
across the channel D. 

Let 

n = number of submarines in the barrier 
p-j^ = P(kill a transiter under I) 

P2 = P(kill a transiter under II) 
qi = 1 - Pi 

q2 = 1 - P 2 

Now, cell width for any single line barrier is related 
to n and D by 




so 

w m 

l-P]^-qi-e =e 



2H 




Alternative I 
The Stacked Barrier 




Alternative II 
The Packed Barrier 



Figure 3 

The Packed and Stacked Barriers 
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and assuming independence among stacked barriers 



1 






where there are k stacked barriers, each with n^ submarines 
(i = 1, . . . ,k) . Now 




and 



= P. 



Thus there is nothing probabilistic to distinguish between 
the two barriers under independence. 

Assuming,' on the other hand, that there exists some 
dependence among the k stacked barriers, then either type 
of dependence is possible. The type of dependence results 
from the nature of the opposing forces. There are only 
four possible results from a transit across the barrier. 
They are the following states : 

1 Barrier killed transiter survives 

2 Transiter killed barrier survives 

3 Both submarines killed 

^ No submarines killed 

In the case of 2 and 3, the type of dependence between the 
barrier in question and subsequent barriers is unimportant 
since the transiter submarine has been killed, and 
thus will not encounter any subsequent barrier. 
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In the case of 1, it is likely that the combat between the 
two submarines would alert the next barrier, or possibly 
several barriers, and allow for maximum readiness and 
optimum placing. In the case of 4, whatever special 
circumstances allowed the transiter to evade the first 
barrier would be likely to carry over into the second 
barrier. These arguments are expressed analytically as 
follows. Let 

= P(enter state jjleave state i) 

1 = 1 , 2 , 3, 4 
5 = 1 , 2 , 3, 4 

where one state transition occurs upon transiting one 
barrier. Then in these cases 

^12 ^ ^Ij ^ ^ 1,3,4 

and 

^ J “ 1,2,3 

Remembering that states 1 and 4 correspond to transiter 
not killed, and states 2 and 3 correspond to transiter 
killed it can be seen that if the first encounter results 
in state 4, then the transiter will be most likely to 
avoid being killed by a state 4 result in subsequent 
barriers. If the first encounter results in state 1, then 
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I 



i 




the translter Is more likely to be killed by a state 2 
result In subsequent barriers. Thus, depending upon the 
likely nature of the first encounter, the decision can be 
made whether to stack or pack barriers. 
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III. PARAMETRIC PREDICTION OP SWEEP WIDTH 



A. DESCRIPTION OP INPUTS 

The applicability of the Random Search Model being 
established, the next step is to find a method for predicting 
the value of W in Eq. tl). Since the analysis is intended 
to serve as an aid to decision making, it is desirable to 
use variables in the prediction whose relationship to W is 
simple, and which are under the control of the decision- 
maker, or tactical commander. 

In submarine warfare, an important measure of effective- 
ness of a submarine is its acoustic quietness, called 
PIGURE-OP-MERIT (POM). POM measures the self noise of the 
submarine, as well as the signal level of the opponent. 

It combines these quantities into one, expressing the 
quietness of the submarine relative to its opponent. 

The POM used in this paper is as follows: 

POM = - DT 

S 0 



Lg = 10 log^^ClO 




+ 10 



L ' /lO 



) 



Lg = 10 log^pClO 



, ea 



+ 10 




+ 10 




) 



L = Radiated Noise 
s 

= radiated machinery noise 
sm 

L = radiated velocity noise 

S V 
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Self Noise 



= ambient self noise 
ea 

** machinery self noise 

L „ = velocity self noise 
ev 

DT = Detection Threshold 

Velocities for computation were taken as Vg, the 
barrier patrol speed, and V^, the transit speed. After 
detection, velocities changed for both submarines, but no 
account of this was taken, since this represented mainly 
kinematics. 

B. THE ESTIMATING RELATIONSHIP 

The dependent variables to be predicted were and Wg 
(sweep width of transiter and barrier, respectively). 
Independent variables were assumed to be POMg and FOM,p 
(FIGURE OF MERIT Barrier and Transiter, respectively). A 
linear relationship was found to be the best predictor of 
Kill Sweep Width, with form: 

= a + b FOM^ + d FOMg a,b,d,f,g,h > 0 

( 2 ) 

Wg = f - g FOM^ + h FOMg 

Both kill sweep widths were tested, that for the 
barrier (Wg) and the transiter (W^), using stepwise linear 
regression. In each case, it was found that the most 
significant variable for predicting sweep width was the 
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FOM of the submarine to whom the W belonged. That is, the 
most significant variable in predicting was POM^. Both 
FOM's were found to be significant for both sweep widths. 

The coefficients preceding the variables in Eq. (2) 
were all as might reasonably be expected, except for d, 
the coefficient of FOMg. It would seem that the sign of d 
should be negative, so that Improvement in the performance 
of the barrier submarine would degrade the performance of 
the transiter . This would ordinarily be the case, except 
for the implications of the tactical assumptions. Since 
the transiter does not seek combat , all opportunity for 
such combat arises by way of the barrier detecting the 
transiter. Thus, as the barrier improves in performance, 
more opportunity to attack arises for the transiter. 

In order to ensure adequacy of the regression, two 
additional tests were conducted. For each, the difference 
between the predicted values of sweep width and the actual 
values, called residuals, were found. These residuals were 
tested using Analysis of Variance. 

The first test concerned the hull types. The question 
tested was whether there was any information associated 
with the Hull/Sonar combinations which was not conveyed by 
means of FIGURE OP MERIT. The results indicated that there 
was no reason to believe that any additional Information 
about the submarine was likely to be of any use, beyond 
FIGURE OF MERIT. 
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The second test was whether the patrol speed (Vg) or 
the transit speed (V,p) had any effect beyond that conveyed 
through POM. The results indicated that the patrol speed 
of the barrier submarine was not significant, whereas the 
transit speed (V^) was significant. This Is due to two 
factors. First, the faster the translter travels, the less 
time he is exposed to detection. Second, the translter is 
less susceptible to intercept at higher speed. This 
remaining effect due to is a drawback in the model, 
since the best possible prediction is not being carried 
out. In view of the excellent fit of the linear regression 
(see the Statistical Appendix) the model seems adequate 
for most needs. It is possible that the additional compllca- 
aion of the model or the increased need for model inputs 
would make a more refined model undesirable. 

C. A SECONDARY EFFECT OF THE ANALYSIS 

In the course of the statistical analysis of the 
simulation data, an important result was obtained, above 
and beyond the requirements of data compactness and 
manageability. The result was the detection, from 
statistical considerations, of programming errors in the 
simulation. Tactical simulations are always large, in a 
computer code sense, and difficult to completely "de-bug". 

In this instance, a logic control variable had been 
included, to allow changing of torpedo P(kill) during a 
run. The program was set up initially using dummy P(kill)’s 
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for all torpedoes. When the final data runs were made, 
the logic control variable was set to inject the correct 
P(klll)'s for barrier types 1 through 4, but was erroneously 
omitted for barrier types 5 through 7. The omission was 
not noticed until the statistical analysis began. As a part 
of the statistical analysis, the residuals (Actual W minus 
predicted W) were plotted for all submarine types. It was 
noticed that barrier types 1 through 4 were all performing 
poorer than predicted, and barrier types 5 through 7 were 
performing better than predicted. This discontinuity was 
not due simply to different Hull/Sonar types, because one 
Hull type and one Sonar were found In both groups, although 
not together. When the programmer saw this effect he 
correctly deduced the error and re-ran the entire set of 
barrier types 5 through 7. The resulting residuals were 
again plotted, and no longer displayed the discontinuity. 

It is unlikely that the error would ever have been discovered 
without the statistical analysis. 
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IV. CONCLUSIONS 



The philosophy pursued herein, that a model is valid 
if it is statistically supported is a profitable philosophy 
only insofar as it aids analysis. If a model is reasonable 
from an analytical point of view and also supportable 
statistically, then the best solution has been arrived at. 

It is the case here, as in many Instances, that the model 
chosen was not uniquely applicable, nor were its requisite 
assumptions completely satisfied. Hopefully, any doubts 
about applicability were dispelled by the compelling 
statistical evidence. 

The method of inversion of results to obtain unknown 
inputs, used here to determine W, is also somewhat 
unorthodox. The resultant estimates of Inputs are naturally 
only as good as the original results. Assuming an acceptable 
level of validity, these Inverted quantities may be compared 
to engineering or similar estimates, probably to the benefit 
of both the engineering and tactical analyst. Perfectly 
valid and useful models are no doubt frequently rejected 
because the Improper estimation of input quantities causes 
inexplicable differences between actual and predicted 
results. It seems reasonable, however, that statistically 
verifiable consistency in quantities Inverted by way of a 
model would argue in favor of the model, and cause 
re-examination of the engineering estimate in case of 
disagreement. Conversely, models which are supposed to be 
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invalid due to assumptions being clearly violated but which 
are found to be statistically supportable can often be 
extended to include the case in question by further 
considerations of the problem. The phenomenon of clipping 
developed here is an example of this potential. 
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APPENDIX A 



STATISTICAL APPENDIX 



1. Validation of the Model 
The formula 



_ W 

p = P(kill) = 1 - e ^ 



W ® kill sweep width 
C = cell width 



was inverted to yield 



W = - C In (1 - p) 



which then becomes a maximum likelihood estimator of W. 
The hypotheses were 

W 




Hq is false 



or equivalently 

Hq: W=-cln(l-p) 

^1' ^0 false 
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The test was by analysis of variance using a Completely 
Randomized Factorial design with treatments of barrier type, 
translter type and cell width at seven, four and three 

levels, respectively. Submarine types were taken as 
fixed factors, cell width as a random factor. Two separate 
tests were conducted, one for Wg and one for (barrier 
and translter kill sweep widths, respectively). P ratios 
and significances were as tabulated In Tables I and II. 

Based on the significance of the cell width treatment, the 
null hypothesis was accepted at significance level of .05. 

2. Parametric Prediction of Sweep Width 
The estimating relationships 

= a + b FOM^ + d FOMg 
Wg = f - g POM^ + h FOMg 

were tested by linear regression. The test statistics for 
were : 

F = 8.9 
= .787 
n = l68 

The regression was thus significant at the level of .05. 

The t statistics associated with the constants b and d 

\ 

were l8.71^ and 2.98 respectively, both significant at .05. 
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TABLE I. 


ANALYSIS OF VARIANCE 


OF CELL WIDTH 


FOR W^ 


Treatment 


Degrees of 


F 


Significance 




Freedom 


Ratio 




Barrier 


6 


80.1 


a 


Transiter 


3 


153.9 


a 


C 


2 


2.015 


b 



TABLE II. 


ANALYSIS OF 


VARIANCE 


OF CELL WIDTH 


FOR Wg 


Treatment 


Degrees 


Of 


F 


Significance 




Freedom 




Ratio 




Barrier 


6 




76.9^ 


a 


Transiter 


3 




181.621 


a 


C 


2 




1.6515 


b 



Notes : 

a « .001 

b > .25 
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In each regression, data Included each of twenty-eight 
submarine pairings at two separate values of Vg and three 
separate values of making a total of l68 data points. 

In each regression there were no other variables ever 
tested as Independent variables, thus the significance 
level reported is valid. 

Upon completion of each regression, the residuals were 
tested for any variability remaining which could be removed. 
This was done by considering two analyses of variance of 
each set of residuals. 

The first post-regression analysis of variance used 
barrier and translter hull/sonar type as treatments at 
7 and 4 levels respectively. The hypotheses were 

Hq: Hull/Sonar type has no effect on the residual 

H^: Hq false 

Test design was Randomized Block Factorial Design with 
a block consisting of the six distinct velocity pairings. 

The two treatments were considered as fixed. P ratios and 
significances were as tabulated in Tables III and IV. Based 
on the resulting significances, the null hypothesis was 
accepted. The overall significance was .05. This implies 
an individual significance of .013, where 

P(No test value outside critical region) = (1-a)^ 

thus, (l-a)^=.95 

a = 1 - (.95) ^ = .013 
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TABLE III 



ANALYSIS OP VARIANCE OP RESIDUALS OP W, 
(SUBMARINE-TYPE TREATMENTS) 



Treatment 


Degrees of 


P 


• Significance 




Preedom 


Ratio 




Barrier 


6 


:i\8 


.8173 


Translter 


3 


1.04 


.3991 



TABLE IV. 


ANALYSIS OP 


VARIANCE OP RESIDUALS 


OP Wg 




(SUBMARINE-' 


lYPE 


TREATMENTS ) 


Treatment 


Degrees 

Preedom 


of 


P 

Ratio 


Significance 


Barrier 


6 




.6 


.7267 


Translter 


3 




- 4.13 


.0216 



The second post-regression analysis of variance used 
Vb and as treatments at two and three levels respectively. 
The hypotheses were: 

Hq : Vg has no effect on the residuals 

Hq false 

and 

Hq! has no effect on the residuals 

H^: Hq false 

For this experiment, there were four separate tests, 
one for each unit's velocity, and thus two for each 
regression. The likelihood is that one or more of these 
null hypotheses might be true with no implications for the 
others. Experiment design was again Randomized Block 
Factorial with each effect now taken as random. Each block 
consisted of 28 distinct submarine pairings. F ratios and 
significance levels are tabulated in Tables V and VI. 

Based on the resulting significances, the velocity of the 
transiting submarine is significant for the probability of 
its detection and kill by the barrier submarine. This 
result is probably due to the resulting reduced time 
available to the barrier submarine to detect, and to the 
increased difficulty in closing with and killing given 



detection. 



TABLE V. 


ANALYSIS OF VARIANCE OP 
(VELOCITY TREATMENTS) 


RESIDUALS 


OP 


Treatment 


Degrees of 
Freedom 


P 

Ratio 


Significance 




1 


3.5 


.20 


v „ 


2 


7 .^ 


.12 



TABLE VI. 


ANALYSIS OP VARIANCE OF 
(VELOCITY TREATMENTS) 


RESIDUALS 


OP Wg 


Treatment 


Degrees of 


F 


Significance 




Freedom 


Ratio 






1 


6.1 


.13 




2 


24.7 


.039 
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APPENDIX B 



THE CLIPPING PHENOMENON 



Suppose that a searcher has Lateral Range Curve p(x) 
(vanishing after R) the probability of detecting a target 
passing at lateral range x. Define 

R 

W=2/ p(x)dx 
0 

Then if the target is uniformly distributed across width C, 
centered at the observer. 



Now, however, suppose that the curve p(x) vanishes 
outside a range R greater than C, and that the searcher is 
constrained to remain within width C and thus may not 
pursue nor take any action with respect to'-a target outside 
of C. Let the searcher's position from the edge of the 
channel be denoted u (Figure 4). 

In this case, the probability of detecting a target 
with lateral range (closest point of approach) distributed 
uniformly across C, for a searcher at u is 



P(detect) = ^ 



R < C 



P(det) = / p(x) ^ dx + 
0 ^ 



0 



/ p(x) ^ dx 




No clipping, searcher centered in channel 




Clipping, searcher at u 



Figure 4 
Clipping 
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denoted hereafter as 



P(det) = i {P(u) - P(C-u)} = i Wg(u) < i W 



where 

P(u) 



u 

/ p(x) dx 
0 



Now, If the searcher Is In turn uniformly distributed 
in position u across C, the probability of detection is 

1 ^ 1 ■ 

P(det) = 7=r / P(x) - P(c-x) dx “ tt W 
0 0 b a 

Now, if the searcher is located in the center of the 
channel. 




and so 

P(det) = ^ (|) 

Clearly, if p(x) is decreasing In x for x positive, and 
symmetric then 




We(§) 




Thus against a uniform target, probability of detection 
Is greatest for a central position, and Is diminished by 
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1 



patrolling. Of course. If patrol Is not conducted, the 
possibility arises that barrier station location may become 
compromised, so a patrol is clearly necessary. 

Let p(x) be of the form 

p(x) = e”^^ a > 0 



which is reasonable for passive sonar detection (see Ref. 1), 
but intended only as an illustration. If the searcher is 
at u = ^ (centered in the channel) 

P(det) = §• Wg (|) = I (1 - e ^ ) 

If the searcher is uniformly distributed across the channel 

2C 

P(aet) = 1 = « (1 - ^(1 - e' «) ) 



These quantities are graphed in Pig. 5 as a function of 
W 

along with the quantities 

W 



and the Definite Range Law (Ref. 1) 



P 





1 

1 
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l.Or ^ (Detection) / Bandom Search Law 
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The agreement Is seen to be excellent among the first three in 
areas of mild clipping, and when there is excessive clipping, 
the Random Search Law (Equation (1)) is somewhat optimistic. 
This contrasts sharply with the obvious Inaccuracy of the 
last quantity, the Definite Range Law. It should now be 
fairly clear- that the amount of clipping depends a great 
deal on the shape of p(x), the lateral range curve. 

Lacking any better information, the slight degradation 
implied in the Random Search Law is quite probably suffi- 
cient to a reasonable degree of accuracy. If, as in this 
analysis, the data support the model, then there should be 
no doubt as to its reasonableness. 

Most cases of barrier patrol will probably be found to 
lie between the cases of the searcher centered in' the 
channel and the searcher distributed uniformly across a 
cell. The problem treated here was within the domain of 
excellent agreement (moderate clipping), and hence the 
law should reasonably hold, since the Lateral Range Curve 
should have been roughly exponential in x. Thus the 
applicability of the formula 

- ^ 

P(klll) = 1 - e ^ 



is established. 
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